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Abstract 

A new sandpile model is studied in which bonds of the system are inhibited 

for activity after a certain number of transmission of grains. This condition 

impels an unstable sand column to distribute grains only to those neighbours 

which have toppled less than m times. In this non-Abelian model grains 

effectively move faster than the ordinary diffusion (super-diffusion). A novel 

system size dependent cross-over from Abelian sandpile behaviour to a new 

critical behaviour is observed for all values of the parameter m. 
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The concept of Self- Organized Criticality (SOC) was introduced to describe how a sys- 
tem, starting from an arbitrary initial condition may evolve to a scale free critical state 
following some specific dynamical rules while under the action of repeated external per- 
turbations [1]. Naturally occurring physical phenomena like sandpiles [2], forest fires [3], 
river networks [4], earthquakes [5] etc. are argued as systems showing SOC. To demonstrate 
the idea of SOC a simple model known as the 'sandpile' model was introduced in which a 
stochastically driven cellular automata evolves under a non-linear, diffusive, self-organizing 
mechanism leading to a non-equilibrium critical state [1]. 

At present many different versions of the sandpile model are available. However precise 
classification of various models in different universality classes in terms of their critical 
exponents is not yet fully complete and still attracts much attention [6]. Among the different 
models most widely studied is the Abelian Sandpile Model (ASM) in which many analytical 
[7] as well as numerical [8] results are known. Some efforts have also been given towards 
the analytical calculation of avalanche size exponents [9,10]. Secondly a Two-state sandpile 
model with stochastic evolution rules was also studied [11] which was initially thought to 
belong to the same universality class as that of ASM [11,12] but later claimed to be different 
[13]. 

We consider a situation in which an intermediate time scale is associated with every bond 
of the system. Each bond allows only a certain number of grains to cross from its one end 
to the other and after that it has a dead time and cannot support any further traffic untill 
a new avalanche starts. This dead time is much greater than the time scale of avalanche 
propagations but much less than the input rate of grains. We call this model as the 'Sandpile 
Model with Activity Inhibition' (SMAI). 

Similar to different sandpile models we also define our model on a regular lattice with 
open boundary. Non-negative integer numbers (hi) assigned at the lattice sites represent the 
heights of the sand columns. Sand grains are added at randomly chosen sites by increasing 
the h values by unity : hi — > hi + 1. The possibility of a sand column becoming unstable 
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arises only when the height hi becomes greater than a threshold value h c . Such a column 
becomes unstable only if the number rii of nearest neighbour sites which have toppled less 
than a pre-assigned cut-off number m within the same avalanche is found to be non zero. An 
unstable column immediately topples and distributes one grain each to all the rii neighbours 
: hj — > hj + 1 (j = 1 to rii). The sand column decreases by the same amount : hi — > hi — rii. If 
rii = 0, the sand column does not topple and its height though greater than h c is considered 
stable. In an avalanche sites can topple a maximum of m times. This implies that in the 
limit of m — > oo our model converges to ASM. Recently a Stochastic Sandpile Model (SSM) 
has been studied in which sand columns having heights greater than the threshold are also 
considered stable [14]. 

One unit of time within an avalanche consists of the following intermediate steps : (i) a 
list of all sites where hi > h c is made (ii) values are calculated for each site i (iii) all sites 
with non-zero rii values are toppled in parallel. 

We first consider the case where the cut-off in the toppling number m — 1. Here the 
toppling front moves outwards and grains always jump only in the outward direction and 
do not fall back. Therefore compared to the random walk analogy for the movement of the 
grains in ASM [15] in our model grains move faster than diffusing particles. This is indeed 
reflected in the average cluster size < s >~ L 13 " where f3 s = 1.62 (reported below). This 
implies that the displacements 1Z of the grains in our model grows with time T as 1Z ~ T u 
with v = 1/1.62 = 0.62 which is faster than diffusion (super-diffusion). 

Zhang had studied a scaling theory of the sandpile model in which the toppling front 
grows as a (d — 1) dimensional surface in the d dimension [16]. In this analysis, where 
multiple topplings were ignored, the exponent for the avalanche size distribution was possible 
to calculate analytically in all dimensions. Since in SMAI, a single toppling front moves 
outward and multiple topplings are forbidden for m = 1, we expect that SMAI may be a 
correct realization of the Zhang's theory [16]. 

Unlike ASM our model turns out to be non-Abelian. Different steady state configurations 
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are obtained on dropping grains at same locations but following different sequences. We pick 
up a stable configuration in a 2 x 2 cell and drop one grain each at the two opposite corners 
one after the other (figure f). On reversing the order of dropping a different stable state 
configuration is obtained. We notice that non-abelianity is effective only when avalanche 
cluster sizes are greater than one. 

However, it is easy to see that just like in ASM, we cannot have the forbidden subcofig- 
urations (FSC) anywhere in the lattice. An FSC is defined as the subset of connected sites 
for which at each site the height is less than its coordination number in the subset [7]. In 
SMAI also two neighbouring sites whose heights are both zero (0 — 0) will never occur in 
the steady state because of the same reason as in ASM that if one topples the other site 
will receive one grain. Similarly a height configuration like (0 — I — 0) is also an FSC. In 
fact, all the FSCs defined for ASM are also forbidden here, and a recurrent config must burn 
completely. This is checked by burning many successive steady state configurations where 
the fire starts from the boundary as well as from those sites with heights greater than h c [7]. 
All the configurations are observed to burn fully. No unburnt configuration is found which 
is obtained by adding grains to a steady state that burns. 

To use the rotational symmetry of the system the sandpile is grown with h c = 3 within 
a circular region of radius R = (L — l)/2 placed on a square lattice of size L x L. In the 
steady state starting from the boundary the average height grows quickly radially towards 
the centre following a power law : < h(r) >= A — B(R — r)~ 5 , where r is the radial distance 
measured from the centre. We estimate A = 2.3904, B = 7.81 and 5 = 0.75 for L = f025 
(figure 2). The average height per site is found to depend on L which on plotting with 1/L 
extrapolates to a value 2.3840 in the limit of L — > oo. Similar analysis yields the fraction 
of sites with different column heights are approximately 2.1xl0~ 4 (h — 0), 0.2421 (h = 1), 
0.3059 (h = 2), 0.3404 (h = 3). Beyond h = 3, this fraction decreases approximately 
exponentially as exp(— ah) where a = 1.64 and adds up to a total of 0.1118. 

The size of the avalanche is measured in three different ways : (i) the total number sites 
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s which crosses the threshold h c (either toppled or not, both counted) (ii) the life time of 
the avalanche (t) and (iii) the linear extent or the radius (r) of the avalanche. Since s, t 
and r are the three different measurements of the same random avalanche cluster, they are 
necessarily dependent variables. These quantities are assumed to depend on one another as 

s-t 7is , r~t 7tr , s~r 7rs . (1) 

The three 7 exponents are connected by the relation 

Its = Itrlrs (2) 

To estimate the exponents 7 ts and j tr we measure the avalanche size s and avalanche 
radius r at every time step t during the progress of each avalanche. The total number of 
topplings up to time t gives the intermediate size s where as the size of the smallest square 
which encloses the cluster gives the intermediate radius r. We estimate 7 <s = 1.64 and 
7i r = 0.83. Since the avalanche clusters are quite compact and has only few small holes it is 
justified to assume that j rs = 2. These values very closely satisfy the equation (2). 

We assume the finite size scaling forms for the probability distribution functions as 

Pi')—- n f.(£;), -*-/.(£). nr)~r-"f r {-^). (3) 

Consequently the cumulative probability distribution F(x) = f£ P(x)dx varies as x L ~ Tx . 
However, in the case of r x — 1, the variation should be in the form F{x) = C— \og{x). 

We plot the data of F(s) in two different ways. In figure 3 we plot F(s) vs. s for system 
sizes L = 65, 257 and 1025 using a log — lin scale. Presence of humps in the large s limit is 
visible for bigger system sizes which reflects the effect of the finite system size on power law 
distributions. However in the intermediate region curves are reasonably straight indicating 
that the exponent r s is likely to be 1. We further plot F( y s)s Ts( - L '~ 1 with s on a log — log 
scale and tune t s (L), the effective r s exponent for the system size L, such that the curves 
become horizontal in the intermediate range of s. All three curves collapse nicely when the 
abscissa is scaled as sL~ 1S2 which implies that a s = 1.62. We show in the figure 4 that the 
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t s (L) values very closely fit to a straight line when plotted with L -1 / 4 . It seems that L -1 / 4 
may be the right leading correction to scaling. The fitted straight line when extrapolated 
to L — > oo gives a value of 1.016 for r s . Similar analysis for the life time distribution also 
leads us to conclude that r t = 1.02, a t = 0.98. 

The radius distribution F(r) is calculated in a slightly different way. We estimate the 
probability that a site at a distance r from the centre of mass of the avalanche cluster is 
a part of the cluster. We take into consideration the 'degeneracy' effect that different sites 
could be at exactly equal distances from the centre. In figure 5 we show a scaling plot 
F(r)L - 20 against rL~ - 86 using a log — lin scale for different system sizes. Here we see a 
much better straight part in the intermediate region. We conclude a value of r r ~ 1. 

The distribution functions follow relations like P(s)ds ~ P(i)dt which imply following 
scaling relations 

t s - 1 = 7t s (r t - 1), r r - 1 = 7t r (r t - 1), r s - 1 = 7 rs (r r - 1). (4) 

These equations imply that if one of the exponents r s , r t or r r is equal to one, rest are also 
equal to one, irrespective of the values of the 7 exponents. Our estimates for the different 
r exponents are very much consistent with these equations. We also observe that the value 
of t s 1 agrees very well with the Zhang's result r s = 2(1 — 1/d) for d — 2 [16]. 

We also assume that the average values of s, t and r varies with the system size L as 

<s(L)>~L /3 % <t(L)>~L A , <r(L)>~iA. (5) 

We plot < s(L) > vs. L on log - log scale for L = 33, 65, 129, 257, 513 and 1025. Slopes 
between successive points are plotted with L~ 2 and extrapolated to L — > 00 limit giving 
/3 S = 1.61. Similar analysis gives j3 t = 0.96 and (3 r = 0.82. 

Using the scaling forms in equation (3) we get following scaling relations for (3 exponents 

as 

(3 s = a s (2-r s ), Pt = <T t (2-T t ), (3 r = a r (2-T r ). (6) 
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With our measured values of (3, a, and r these relations are approximately satisfied. We put 
errors of 0.05 to all our measured exponents. 

Next we study the case when cut-off for the toppling number to > 1. The average cluster 
size < s(L) > is plotted with L on the log — log scale in figure 6 for to — 1, 2, 4 and 8. We 
see that all curves are parallel straight lines with slopes approximately 1.61 for large system 
sizes. However for small system sizes all of them bend and become part of the same straight 
line. Then we plot on the same figure the < s(L) > data for ASM. We get a straight line 
with a slope « 2 which almost overlaps with the bend portions of the curves for different 
to values. We explain this by noting that for every to value our model behaves as ASM for 
small system sizes. In small systems the number of avalanches where sites will topple more 
than to times are very few. However for bigger system sizes the cut-off to will have more 
prominent effects. Therefore for each to values there should be one particular system size 
where the cross-over takes place from ASM to Non-abelian behaviour. The cross-over size 
L c is seen to be roughly proportional to the value of m. We expect that for any to value if 
one works in systems larger than the cross-over size one should get the same set of exponents 
as those in the case of to — 1. 

To summarize, we studied here a new sandpile model where bonds of the system relax 
after a certain number of transmission of grains. This limits a site to topple a maximum of 
to times within the same avalanche. Based on the results of detailed numerical studies using 
improved algorithms we claim a cross-over from ASM behaviour to a new critical behaviour 
at a particular size of the system whose magnitude depends on the value of to. 

We acknowledge with thanks D. Dhar for many useful discussions and suggestions. 



7 



FIGURES 

FIG. 1. Non-abelian property of the sandpile model is shown on a 2 x 2 cell. On the same initial 
stable configuration two grains are added at two different sites but in different orders. Different 
final stable configurations are obtained. 

FIG. 2. The average height profile of the sandpile in a circular region plotted with the distance 
from the circumference of the circle is shown. The slope of the curve is 0.75 and A = 2.3904 is 
found. 

FIG. 3. Log — lin plot of the cumulative probability distribution F(s) for the three system sizes 
L = 65,257 and 1025 (from left to right). The straight portions of the curves in the intermediate 
regions indicates that t s is likely to be equal to one. 

FIG. 4. Plot of t s {L) for different system sizes L = 33,65, 129,257,513 and 1025 with L" 1 / 4 . 
A direct straight line fit gives r s = 1.016 in the L — > oo limit. 

FIG. 5. Scaling plot of the cumulative radial distribution function F(r). Plot of F(r)L - 20 vs. 
r/L 0,86 shows the data collapse for the system sizes L = 65, 257 and 1025. 

FIG. 6. Plot of < s(L) > versus L for m = 1, 2, 4 and 8 of SMAI (solid lines) and for ASM 
(dot dashed line). For each value of m there is a threshold system size L at which the cross-over 
from ASM behaviour to SMAI takes place. 
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Figure 1 : Manna and Giri 
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Figure 4 : Manna and Giri 



4.0 
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